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The possibility that particle production in high-energy collisions is a result of two asymmetric 
hydrodynamic flows is investigated, using the Khalatnikov form of the 1+1-dimensional approxi- 
mation of hydrodynamic equations. The general solution is discussed and applied to the physically 
appealing "generalized in-out cascade" where the space-time and energy-momentum rapidities are 
, equal at initial temperature but boost-invariance is not imposed. It is demonstrated that the two- 

(N . bump structure of the entropy density, characteristic of the asymmetric input, changes easily into 

^ ' a single broad maximum compatible with data on particle production in symmetric processes. A 

, possible microscopic QCD interpretation of asymmetric hydrodynamics is proposed. 

Q 

(N 

I. INTRODUCTION 

Oh; 

^ . The standard phenomenological description of heavy-ion reactions is based on the hydrodynamic evolution of the 
system, starting from rather early times. This approach allows to explain many aspects of data, e.g. the transverse 
momentum spectra and elliptic flow [1, 2]. 

The hydro calculations mostly concentrate at the central rapidity region and therefore they often leave aside the 
^ ' problem of longitudinal dynamics. The usual approach for symmetric collisions {e.g. for gold-gold collisions) is to 
consider a rapidity symmetric medium which, after some pre-equilbirium stage, evolves following the hydrodynamic 
equations till the freeze-out temperature. The investigated solutions of the hydrodynamic equations are naturally 
. symmetric since they are obtained with initial conditions which are themselves rapidity-symmetric. 
• ' In the present paper we discuss a different hydrodynamic approach where, starting from asymmetric initial con- 
^ ^ , ditions, one considers asymmetric solutions of the hydrodynamic equations. The symmetry of the overall solution, 
'■ in particular of the rapidity distribution of particles, is thus restored by the superposition of two conjugate asym- 
, metric solutions. Such a situation has been advocated to explain (i) the forward-backward multiplicity correlations 
observed in nucleon-nucleon and nucleus-nucleus reactions [3-5] as well as (ii) the sign and magnitude of the directed 
flow observed in nucleus-nucleus collisions [6]. It is also natural in some models [7-10] describing the multiplicity 
distributions in p-p and gold-gold reactions. Finally, it is obviously necessary for an hydrodynamic description of the 
5— I ' asymmetric reactions, e.g. deuteron-gold collision [11]. 

For this scenario to be realized, one has to show that the hydrodynamic evolution is compatible with the existence 
of asymmetric components and with their possible combination to form a reasonable overall symmetric distribution. 
It is the goal of our paper to show that this is indeed the case. We use the simplified but informative framework 
of the (1-f l)-diniensional approximation of the corresponding hydrodynamic fiow. In this context we i) discuss 
the construction of asymmetric solutions, ii) show that the linear superposition of asymmetric components of the 
multiplicity distribution in rapidity is possible to construct in an hydrodynamic framework, more precisely of its 
(l+l)-dimensional approximation by the longitudinal entropy density as a function of rapidity. 

It was shown long time ago by Khalatnikov [12] that the hydrodynamic equations in (1+1) dimensions, considered 
first in the context of high-energy collisions by Landau [13, 14], can be linearized. This is achieved by passing 
from kinematic variables (position and time) into dynamic ones (temperature and rapidity). Recently this idea was 
developped in a modern context [15, 16] and used in the discussion of the longitudinal evolution of the symmetric 
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solutions. The argument of the present paper extensively exploits these new developments and in particular the 
linearity property, which allows to consider the superposition of various solutions^. 

We shall discuss how one can obtain explicit solutions of the hydrodynamic equations which can be in phenomeno- 
logical agreement with the two-source description in terms of the wounded constituent models [7, 8, 11]. In doing so, 
we will more generally describe the full exact solution of the (l+l)-dimensional approximation of a perfect fluid flow. 
Finally, we shall speculate on the possible field-theoretical interpretation of our solutions in terms of early asymmetric 
thermalization in heavy-ion high-energy collisions. 

The plan of our paper is the following: in section II we recall the derivation of the Khalatnikov equation. It allows 
for a general solution of the (l+l)-dimensional approximation of the perfect fluid flow which we derive in section III. 
From that general solution we discuss the existence of conservation laws as shown in IV, and then the definition of the 
initial conditions which are formulated in V. In section VI, we derive and discuss in detail the physically appealing 
family of solutions corresponding to the generalization of the "In-Out cascade " scenario [17], where the space-time 
and momentum rapidities are taken equal at initial temperature. In the final section VII, we give our conclusions and 
discuss the possible microscopic QCD interpretation of asymmetric hydrodynamics. 



II. THE KHALATNIKOV EQUATION 

It is known [12, 14] that one can formulate the problem of (1+1) hydrodynamic evolution with a linear equation for 
a suitably defined potential. In this section we briefly recall these results, recasting the calculations in the light-cone 
kinematic variables, as used in [15, 16] 

z- = z^±z^^t±z = re-^ ^ ^^^±^)=2^i^2d,), (1) 

where r = V z+z~ is the proper time and = | ln(z+/2~) is the space-time rapidity of the fluid. We also introduce 
the dynamical variables: y, the energy-momentum rapidity variable and 6, the logarithm of the inverse temperature, 
deflned as 

y = log u+ = - log u- ; 6 = log(To/T) , (2) 

where = log(u'^ ± v}) are the light-cone components of the fluid velocity and Tq some given fixed temperature. We 
shall choose it to be the initial temperature of the plasma. Since plasma cools during the evolution we have always 

e>Q. 

Using thermodynamic relations, one can recombine the two relativistic hydrodynamic equations coming from the 
conservation of energy- momentum in (1+1) dimension, into two equivalent ones having a direct physical interpretation 
[12]. 

(i) The relation 

d+ {e-'+y) - a_ {e-'~y) . (3) 

implies the existence of a hydrodynamic potential $(^"'", z~) such that: 

d^^{z+, Z-) = To e-^^y = u^T, (4) 

meaning that the flow derives from a kinematic potential. 

(ii) The second independent equation 

d+ (u+s) + d- {u-s) = (5) 

corresponds to the conservation of entropy along the flow. 

In order to combine usefully eqs.(4) and (5), one introduces the Khalatnikov potential [12] 

x{0,y) = ^z+,z-)-z-u+T-z+u-T , (6) 



A recent review of these and other investigations on Khalatnikov eq. can be found in [2]. 
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where arc now considered as functions of {0,y). This is cahed the hodograph transformation expressing the 
hydrodynamic equations as a function of the dynamical variables {0,y) via the Legendre transform (6). Using the 
standard differential relations of a Legendre transform, the kinematic variables are recovered from the Khalatnikov 
potential by the equations 

zHO,y) = ^e'^y {dex±dyx). (7) 

The Khalatnikov potential has the remarkable property [12, 14-16] to verify a linear partial differential equation 
which takes the form^ 

dhiO, 2/) - [1 - c^] dexiO, y) - d^x{0, 2/) = . (8) 

In (8), Cs (denoted also for further convenience) is the speed of sound in the fluid, which will be considered as 

a constant in the present study. 

The Khalatnikov equation has been originally derived [12] for the potential (6) and with applications to specific 
problems. But its range of applicability appears to be much wider. Indeed, the transformation of a nonlinear problem 
in terms of the kinematic variables into a linear one in terms of dynamic variables has tremendous advantages, as we 
shall see further. It applies to other quantities than the Khalatnikov potential; It allows to obtain new solutions by 
arbitrary linear combinations of known ones. Furthermore, primitive integrals and derivatives of solutions are also 
solutions, as shown^ in [16]. 

In order to illustrate the powerfulness of this method, let us observe that if xi^^u) is solution of (8), then the 
potential deflned through (4), but now expressed in terms of the hydrodynamic variables through (7), reads 

^e,y) = ^{z+{e,y),z-{e,y)} = x{0,y)+dex{0,y), (9) 

and thus it verifies also the Khalatnikov equation (8). As a direct consequence, the physical entropy flow as a function 
of rapidity also verifies (8). Indeed, one has [15] (see also [19]) 

^{e,y) = 2soTl e-(s-i)^ a,$(^,y) , (10) 

where we use the thermodynamic relation s = sqT^ = soTqC"^^ for the overall, temperature-dependent, entropy 
density of a perfect fluid, recalling that by definition = 1/g. Since the entropy distribution gives approximately the 
predicted distribution of produced particles, eq.(lO) provides a rather direct link between solutions of the Khalatnikov 
equation and observables. Thus (10) plays a central role in the present study. 

Before going further, it proves useful to use appropriately rescaled variables allowing to put the Khalatnikov equation 
(8) in a simple and symmetric form. Introducing 

t=i^e x='—ly, (11) 

redefining x{t, x) as a function of these reduced variables, and introducing the auxiliary function 

Z{t,x)=e-'x{t,x) , (12) 
the Khalatnikov equation (8) takes the simple form 

{d^-dl)Z = Z (13) 

III. THE SOLUTION 

Equation (13) is nothing else than the 1+1-dimensional Klein-Gordon equation (called also the Telegraphic equation 
[20]). The solution, for arbitrary (Cauchy) initial conditions at t = 0, is [21] 

Z=^{f{x + t)+f{x-t)} + - dsh{s)h{r) + - dsf{s)^ (14) 



^ Note that in this relation there is a sign difference comparing to that of [15, 18], due to the sign-difference in the (^-definition that we 
use in this work. 

^ The paper ref. [16] is based on the family of (1-1-1) dimensional "harmonic flow" exact solutions found in [18]. 
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where 



r = ^/t"^ — {s — xY 

and Io{r) is the modified Bessel funetion of the first kind. 

According to Eqns. (9,10) and using (12) we can express the entropy distribution as 

- = ^(5- l)e-* {{9 + 1)Z + 2gd,Z + {g - l)d^Z} . 



dy 

One also has for the kinematic flow variables, using (7), 
g-l i£H 



2^To 



z + dtZ + cAz 

Z + dtZ- Csd^Z 



(15) 



(16) 



(17) 



IV. CONSERVATION LAWS 



To discuss the conservation laws we observe that it is possible to rewrite the general solution (14) in the form of 
an inverse Laplace transform, implying the identity [22] 



@[t-\s-x\]h yt^-{s-xY ] = ^ 



c+zoo 



(Z7 



7t-|s-x|y'72-l 



(18) 



2m 

with Jie c > 1. Inserting this expression into (14), reversing the order of integrations and using the derivative relation 



^ {e[t-|s-a;|]/o \yjf^-{s-xf\ } ^ e[t-|s-x|]|^/o \sj¥^^^\ + 8\t-\s 



d 



one sees that the solution (14) can be written in the form 



c+icxD 



-7* 



1 



+00 



ds [/i(s) +7/(s)] e-l'*-^l\/^. 



It is not too difficult to show that (18) and (19) do indeed satisfy (13). 
Starting from (20) one can derive a sum rule: 



/ + OO !■ 
dxZ{t, x) = e"* / 
-00 J — 



J — ( 



I dx 

-00 



h{x) - fix) 



where the integral 



POO pc+ioo 1 -ft c+oc 

/ dxZ{t,x)= TT--^^ dx[h{x)+^f{x)] 



was evaluated using the residue theorem. 

From (21) one can obtain two conservation laws: 



/+cx> r+oo 
dx[Z{t,x)+dtZ{t,x)]= dx 
-00 J —00 



h{x) + fix) 



(19) 



(20) 



(21) 



(22) 



(23) 



/+00 
-00 



dx [Zit,x) - dtZit,x) 



"^dx 



(24) 



Comparing (23) with (16) one can see that (23) implies that the total entropy is independent of i, i.e. independent of 
the ratio Tq/T. This means that the total entropy does not change during the evolution of the system. Moreover, since 
the freeze-out temperature may safely be considered as constant, it also means that the total entropy is proportional 
to Tq , as expected from rules of thermodynamics. 

Interestingly enough, equation (24) furnishes another independent conservation law for the 1 + 1-dimensional hy- 
drodynamics of a perfect fluid. Probably due to the (l-l-l)-dimensional origin of the Khalatnikov equation, it would 
be interesting to know its physical interpretation. We checked the validity of both conservation laws in our numerical 
applications. 
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V. INITIAL CONDITIONS 



To discuss the consequences of the Khalatnikov eq., one has to determine, in terms of physical quantities, (i) initial 
conditions at i = (i.e. the functions f{x) and h{x)) and (ii) the final value of t at which the solution is to be 
evaluated. To this end we observe that at t = we obtain from (14) 

Z{x) = fix); dtZ{x) = h{x); d^Z = f"{x) + f{x) . (25) 

Introducing these formulae into (16) and (17) we have 

,_o ^ {-^^^^ ^ ^^""^ ^ %^^"^^'^} ' ^^^^ 

and 

roix) = ^Viim + KxW - [f'ixW/gy, Voix) = N^x{=y) + hog{ + ^f] + "^^^g j (27) 
2To ff-1 2 lf{x) + h{x)-Csj'{x)} 

where {To{x),r]o{x)) = {t,-!])]^^^. 

To determine the final value of t we have to fix the physical conditions at which the system changes into free-moving 
hadrons. It seems reasonable to assume that this happens at a fixed freeze-out temperature T = Tf. Then the final 
value of t at which the physical solution should be evaluated is 

t^i/=log^+log^ . (28) 

We thus conclude that the initial temperature To is needed to be defined in order to determine the final value of t at 
which the physical solution is to be evaluated. 

If To is fixed then eq.(27) determines the proper time tq{x) at which the evolution starts and eq.(17) gives the 
proper time at which freeze-out takes place. To illustrate these features, in the next section we discuss a simple while 
physically appealing case, when f{x) = 0. 



dS 
dy 



VI. GENERALIZED IN-OUT CASCADE 



When f{x) = 0, we have 

Z=l r^'ds h{s) loi^/t^ - (s - x)^) . (29) 

^ Jx-t 

with t = tf given by (28). Using this we also derive 



dtZ = l{h{x+t) + h{x-t)} + \ 12^8 his) J' ^)! ) (30) 

and 



d^Z = l{h'ix+t) - h'{x-t)} + \ {h{x^t) + h{x-t)} + \ 12 ds h{s)^^ 
where we have used the formula 



(31) 



dla{r) A th{r) t , 

^ = I'{r)- = — ^ - when r -» . (32) 
dt °^ \ r 2 ^ ^ 

To see the physical meaning, we observe that when (29) is introduced into (7) one has 

z^{0,x) = !^^ h{x)^ u^; Vo = y\t=o, (33) 

i.e. at f = the position is aligned along the velocity . This can be interpreted as a generalization of the in-out 
cascade [17] for which r] = y at all temperatures leading to boost-invariance hydrodynamics (or "Bjorken flow"). In 
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FIG. 1: dS/dy for 5 growing values of To/Tf = 2'' ,k = Q ^ A. Left: asymmetric one-component model; Right: symmetric 
two-component model. Here one takes c = 0.4 (see text). 




.+ ,+* -tfc % \ 




c = 0.5 




FIG. 2: Two-component dS/dy for 5 growing values of To/Tf = 2*^, fc = — > 4. Left: strongly bumped asymmetric input 
(c = .3); Right : weakly bumped asymmetric input (c = .5). 



our generalized in-out cascade, the relation rj — y is valid at the initial stage t = only. The boost invariance is 
violated (unless h{x) is a constant) because the ratio u/z depends on rapidity. 
Eqs. (26) and (27) reduce to 



dS 
dy 



h{x); To{x) = h{x); r]o = y 



(34) 



To show a specific example, let us consider the following simple parametrization of the intitial entropy distibution 
atT = To(i = 0) : 



h{x) = (1 - x/XmaxT (1 + -f /x^axf exp [(x/x„ 



l)/c] 



(35) 



where Xmax = tmax = logTo/T/ and a,b,c are constants characterizing the initial entropy profile. In the plots we 
keep a = 6 = 2 for simplicity and vary c = {0.3, 0.4, 0.5}. 

In Fig. 1, we display the resuhing curves corresponding to the entropy vs. rapidity distribution for 5 temperature 
values from the initial one to the freeze-out, namely for To/Tf = 2*^, fc = — > 4. It can be seen that the hydrodynamic 
diffusion towards the tails of the one-component distribution (Fig. 1, left) is responsible for the smooth disappearance 
of the two-bump initial structure of the symmetrized distribution (Fig. 1, right) obtained by linear combination. The 
initial bump structure may be easily modified by changing the form of h{x). This also modifies the evolution of the 
entropy distribution as shown in Fig.2. 
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y y 

FIG. 3: Two-component proper-time distribution t{x) for 5 growing values of Tq/T/ = 2*, fc = 4. 
Left: asymmetric one-component model; Right : symmetric two-component model. Here one takes c = 0.4. 

It is also interesting to display some characteristic kinematic features of the flow solutions, both for the asymmetric 

and symmetric cases. Note that for those features, the symmetric configuration is not a linear symmetrized combina- 
tion of both components, due to the nonlinear character of the solution for kinematics (17). The r-distribution can 
be found in Fig. 3. 

Concerning the kinematic properties of the flow one can make the following comments. 

First, as seen in the left panel of Fig. 3, the distribution of proper times is signiflcantly asymmetric. This means 
that, at any given y ^ 0, the two asymmetric components frcczc-out at different proper times. This observation 
supports the idea that the sets of particles emitted by these two components are indepedent of each other (as required 
by the analysis of the long-range correlations [3-5]). 

If the two components arc mixed up during the evolution, the frcczc-out time is very different from the previous 
case. This is illustrated in the right panel of Fig. 3 where one sees that, for a fixed temperature, the proper-time 
of the symmetric solution is close to constant in the central rapidity region, at least at not-too-large values of the 
ratio To/Tf. This is consistent with an approximate Bjorken flow [17]. One also observes that t(j/) increases with 
increasing ratio To/Tf, as expected. Note that the curves for To/Tf = 1 (the lowest ones in the Fig. 3) represent the 
proper time at which the evolution starts. 

We have also evaluated how much the space-time rapidity r] deviates from y along the evolution. As shown in Fig. 
4, we have found that the difference r] — y is small, of the order of at most few %. This conflrms that the the evolution 
of the system we are considering is close to the Bjorken flow. All these features show that the flow is rather smooth 
and laminar. We observe a slight but systematic tendency to obtain rj < y, growing for large space-time rapidities, 
which indicates that the flow is going "outwards" , i.e. moving apart on both sides from a kinematic tangential flow 
which would correspond to the original "In-Out" cascade mechanism. 

VII. SUMMARY AND COMMENTS 

In summary, we have investigated the longitudinal hydrodynamic evolution of a one-dimetional perfect fluid (with 

a constant sound velocity). Using the Khalatnikov equation in the light-cone variables, derived earlier [15, 16], the 
problem could be reduced to the discussion of solutions of the (H-l) Klein-Gordon equation (called also "Telegraphic 
equation" [20]). Using the well-known general form of these solutions, we discussed their consequences in case the 
initial conditions are asymmetric, as necessary e.g. for proton-nucleus or asymmetric nuclei collisions and as also 
suggested by some models [8-11]. 

Wc have studied in some detail only one simple but physically appealing example which can be considered as 
a natural generalization of the Bjorken flow [17] which, as is well-known, give a boost-invariant expansion of the 
fluid. To relax boost-invariance, we assume that although at the beginning of expansion z"^ ~ (c.f. (33)), the 
proportionality coefficient is not a constant (as in the case of Bjorken flow) but it is a function of y. This violates 
boost invariance and thus, among other rapidity configurations, it is possible to study asymmetric ones. The results 
of our study may be synthesized as follows: 

i) Our main conclusion is that the hydrodynamic expansion broadens rather effectively the asymmetric bumpy 
structure of the entropy density (taken for the input). Consequently, for symmetric collisions (as, e.g. p—p or Pb—Pb) 
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FIG. 4: ri—y for 5 growing values ofTo/Tf, marked in the figure; (c = A). Note that r) — y = for To/T/ = 1. Left: asymmetric 
one-component model; Right: symmetric two-component model. 



where the asymmetric input has to be symmetrized, the two-bump structure of the entropy density, characteristic 

of the asymmetric input, changes easily into a single broad maximum. Thus we have shown that even a strongly 
asymmetric input can be compatible with data on particle production in symmetric processes. 

ii) We have also studied the properties of the flow. The most important consequence of the initial conditions 
violating boost invariance is that the proper time at which evolution begins is not constant but depends on rapidity. 
In the simple case we considered, it is even downright proportional to the initial entropy density. Also the freeze-out 
time depends on rapidity, although the; evolution has the tendency to smooth out this dependence. 

iii) The rapidity dependence of tq and r/ is especially interesting for strongly asymmetric initial conditions. In this 
case also the proper time is small in the region of small density and large in the region of large density. Consequently, if 
two asymmetric sources arc present (as is the case for symmetric collisions) they freeze out at different times (except at 
cm. rapidity close to 0). Therefore one may expect that they produce two bunches of particles which are indepenent 
of each other. 

iv) We have also studied the evolution of the difference ry — y (which vanishes at i = in our simple example). It 
turns out that this difference remains small (below of 2-3 % of y) throughout the evolution (except at the very edge of 
phase-space where it may reach 10 %). Thus we conclude that the flow has tendency to be laminar and smooth, not 
very different from the Bjorken picture. It shows some systematic tendency to behave "outwards" i.e. moving apart 
on both sides from a kinematic tangential flow which would correspond to the original "In- Out" cascade mechanism. 

v) As an interesting by-product, we have foimd two general identities valid for the (l+l)-diniensional approximation, 
corresponding to conservation laws along the temperature (and thus proper-time) evolution of the flow. One of them, 
see eq.(23), can be interpreted as the conservation of the full entropy. The other one, see eq.(24), up to our knowledge 
is new, and still asks for a physical interpretation. 

Let us add a comment on the microscopic interpretation of our results for high-energy heavy ion collisions where the 
hydrodynamic set-up proved strikingly pertinent for the description of the Quark-Gluon Plasma phase of QCD [1, 2]. 
It is indeed extraordinary that hydrodynamics describes successfully and almost unambiguously this complicated 
process, belonging to the strongly coupled regime of QCD of which one does not know yet a microscopic description 
(although the string theory gives useful indications through the gauge/gravity duality in a kinematic coUisional 
framework [23], but for yet only supersymmetric gauge fleld theories). 

As is well-known, the application of hydrodynamic description to heavy ion collisions requires a rather rapid 
thermalization of the system of quark and gluons created in the; collision. This is not easy to reconcile with the idea 
that this system is described by perturbative QCD which is the usual scenario considered in this context**. The 
perturbative QCD evolution of the wave-function of a boosted projectile generates a set of soft quantum fluctuations 
(gluons and q — q pairs) which develop as a virtual cascading process. This is revealed, for instance, when these a 
priori virtual fluctuations become observables, as in the case of the "hump-backed plateau" [27] generated by a jet. 
Similarly, such virtual fluctuations play an important role in the theoretical understanding of deep-inelastic reactions. 



* For other approaches, see e.g. [24-26]. 
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The problem arises because this system, considered as the initial state of the incident nuclei is highly anisotropic, 
therefore its thermalization after the collision is difficiilt and would require relatively long time. 

Eventually, however, if the energy is sufficiently large, the soft fluctuations are dense enough in phase-space to 
interact and eventually lead to a saturation mechanism [28] grouping them around some energy-dependent virtuality 
called the saturation scale. Therefore there is a possibility that, due to soft interactions which are long range in time, 
such a system thermalizes (at least approximately) already before the compound system created by the collision takes 
place^. This is a possible model of the initial conditions we considered in the paper. It naturally requires asymmetric 
initial conditions. 

To conclude, our investigation shows that the longitudinal hydrodynamic expansion of the system created in a high- 
energy collision may significantly influence the original rapidity distribution of partons in the colliding projectiles. It 
is therefore important to study this effect in more detail to obtain a reliable information on hadron structure at high 
energies. 
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